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ABSTRACT 


This thesis examines the role of modem artillery as a defensive as 
well as offensive weapon—one which is capable of destroying targets 
and achieving an advantage on the battlefield rather than merely per¬ 
forming traditional tasks of attrition. The thesis demonstrates that the 
new capabilities can be attained. A new concept of artillery weapon 
system—the Trajectory Corrected Artillery Rocket system (TCAR)—is 
analyzed. Results show that this new artillery system, when it contains 
a cluster bomblet warhead, is very effective against infantry in an open 
area, with destruction levels of 50 to 90 percent, but it is not suffi¬ 
ciently effective against fortified and armored targets. However, when 
this system accommodates SFM (“smart") submunitions, it proves 
effective against armored vehicles, with destruction levels of up to 
70 percent. Simulation programs were developed which assessed 
damage levels on a variety of targets. A set of these targets was chosen 
for which a comparison analysis was made between the TCAR and two 
other well-known artillery systems: 155 mm gun and a free-flight 
artillery rocket system. Two parameters were tested: rate of kill and 
marginal cost. Results demonstrate the clear advantages of using TCAR 
over the other systems when accompanied by a cluster warhead. Fur¬ 
thermore, it was found that the TCAR is the currently preferred sys- 


tern suitable for SFM. 
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I. INTRODUCTION 


A. THE ROLES AND CHARACTERISTICS OP MODERN ARTILLERY 

Artillery weapons were never considered main actors on the battle¬ 
field or a destructive force with a imique mission. TTie conventional the¬ 
ory of the role that artillery plays on the battlefield is that it is a weapon 
which supports assault and defense forces—it acts as an attrition and 
suppression weapon assisting main forces. Another important facet of 
the nature of artillery weapons, according to this convention, is that they 
cover area targets (as opposed to pinpoint-accurate antitank weapons). 

However, a modem view is that artillery weapons, in addition to 
their conventional tasks, also act as main forces with unique missions. In 
addition, they have accurate singular destructive capabilities. This new 
concept permits the use of artillery for the following missions: 

• Breaking massive assault forces of superior numbers, and 

• Attacking high-priority long-range targets. 

Modem artilleiy is characterized as: 

• Capable of massive fire power operating in small units 

• Long-range (more than 30 km) 

• Destructive capability for a variety of targets (including “heurd" 

targets) 

• High survivability of launching units on the battlefield 

Newly developed artillery rocket systems were expected to achieve all 
the above requirements, and they are indeed effective. Due to their low 
accuracy, however, th^r were found not to be cost-effective weapons. 
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Newly developed artillery rocket systems were expected to achieve all 
the above requirements, and they are Indeed effective. Due to their low 
accuracy, however, they were found not to be cost-effective weapons. 

B. TRAJECTORY CORRECTED ARTILLERY ROCKET SYSTEM 

The Trajectory Corrected Artillery Rocket system (TCAR) Is a new 
concept of an artillery weapon which Is In a feasibility study stage at the 
Israel Military Industries (IMI) Ltd. TCAR could be defined as a seml- 
guided missile, but we have chosen not to call it such because this family 
of weapons system has an entirely different structure and purpose. The 
system has all the adveintages of the modem artillery rocket systems, as 
mentioned above. By controlling the rocket trajectory, accuracy is sub¬ 
stantially enhanced. As a result, force effectiveness and system cost- 
effectiveness may be significantly Improved. 

1. TCAR Operation Concept 

A rocket equipped with only a simple steering mechanism is 
fired from a launcher. A ground station, located near the launcher, 
tracks the rocket trajectory. At a certain point along the trajectory, the 
ground station computer simulates the expected hit points based on the 
data gathered from this tracking and compares them to the target loca¬ 
tion. A correction command is sent to the rocket, which then performs 
the correction. The tracking system process continues in a closed loop 
until the rocket disappears from line of sight. 

The TCAR Is neither a pinpoint-accurate missile nor an inex¬ 
pensive free-fllght ballistic weapon—it stands between those two cate¬ 
gories. The addition of an Inexpensive (but not perfectly accurate) 
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guidance facility to the free-flight rocket increases the unit cost, but the 
cost is still significantly less than that of a missile. Therefore, the TCAR 
is capable of being employed extensively as an artillery weapon. 

This thesis evaluates the Trajectory Corrected Artllery Rocket 
system’s effectiveness and cost-effectiveness as compared to two other 
artillery systems: 

• Extended range howitzer guns 

• Modem free-flight rocket systems. 

An extended-range howitzer gun is called a Rocket Assisted Projectile 
(RAP) or Base Bleed (BB). Modem free-flight rocket systems are currently 
being used in the US as MLRS and in Israel as IAR160. 

Two types of warheads will be discussed: 

• Cluster warhead with antlpersonnel/antiarmor (ap/am) scatterable 
bomblets 

• “Smart” warheads which contain scatterable “sense and destroy" 
antiarmor/antitank submunitions (This thesis text uses the alterna¬ 
tive name. Sensor Fuzed Munitions (SFM].) 

Models of some typical terrain targets were chosen on which troops, artil¬ 
lery, armored personnel carriers and tanks were located. 

2. More-Detailed Engineering Description Restrictions 

The TCAR concept already exists in the form of an Israeli classi¬ 
fied system, so a complete description of that system is not possible here 
and the engineering details cannot be discussed. However, performance 
data relevant to a generic TCAR system will be used in the models of this 
thesis. 
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C. GOALS OF THIS THESIS 

The goal of this thesis is to generate quantitative data that will be 
used by the system's developer to assess cost-effectiveness in the acqui¬ 
sition of the new system by potential customers (armies). Main points to 
be included are: 

1. The effectiveness which is achieved when artillery rocket system 
accuracy Is enhanced, and 

2. A comparison of the cost-effectiveness of competitive artillery wea¬ 
pons in order to test the feasibility of this new concept. 

It is not within the scope of this thesis to carry out a full combat deploy¬ 
ment analysis of each weapon system. 

D. METHODOLOGY 

The work Includes building simulation programs that will accept, as 
input, all the parameters of those random factors affecting the perfor¬ 
mance of each system. This will result in an evaluation of the effective¬ 
ness of each system on the specific target models. 

The simulation results will determine the number of rounds for each 
system that must to be fired in order to achieve various target damage 
criteria. Cost-effectiveness (i.e., cost per target destruction level) is calcu¬ 
lated accordingly. 

The following additional factors will be considered: 

• Plate of fire 

• Size of the firing units 

• Survivability of the firing units 
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E. INPUT DATA 

The particular systems on which the simulation will be performed 
for comparison with the TCAR are: 

• G155— 155 mm howitzer gun, using as ammunition the RAP (Rocket 
Assisted Projectile) 

• MARS (Multipurpose Artillery Rocket System)—based on the rocket 
system LARI60 (or Its derivatives) that was developed and Is cur¬ 
rently produced by the IMI, Ltd. 

All the data necessary for the simulation was provided by IMI as real 
data. However, because the TCAR is under development, analyses must 
be based on continually updated data for the system as it becomes 
available. 

F. SIMILAR ANALYSIS WHICH HAS BEEN PUBLISHED 

Due to the new military applications of the systems and warheads 
which were analyzed, no comparable studies could be found in the open 
literature. IMI’s Operational Analysis Department and some of the Israel 
Defense Forces operational research branches have done damage analy¬ 
sis for weapon systems that carry cluster bomblet warheads. All those 
analysis results were printed In Hebrew, while some of them, by nature, 
are classified and therefore are unavailable for this thesis. Similar mili¬ 
tary US material (even unclassified) was inaccessible because of the 
author’s status as a foreign student in the US. 

The TCAR as a weapon system is a new concept, so no published 
analysis could be found In the unclassified literature. 

General methods of target damage assessment are described In the 
References 1 through 4. A description of the SADARM (Sense and 
Destroy Submunition) Is given In Reference 5. The SFM, which is 
described elsewhere in this thesis. Is similar. 
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n. THE MODEL AND ITS IMPLEMENTATION 


It was essential to this thesis to write a Monte Carlo simulation 
model In order to evaluate target damage levels for various target types 
using the three alternative systems. Since an area target Is the simplest 
type to be evaluated, the results from the simulation model were checked 
against an exact formulation (using numerical Integration) of the area 
target to determine the accuracy of the simulation model. Once the simu¬ 
lation model generated results comparable to the exact model for the 
area target, the simulation model was modified with several other target 
types and munitions types. Missions on which salvo time duration Is 
considered a deficiency factor and targets that contained scattered 
armored vehicles were also Included. 


A. INTEGRATION METHOD TO ASSESS DAMAGE LEVEL OF 
HOMOGENEOUS TARGET 

1. Damage Assessment Model of a Homogeneous Target 
a. Lethal Area qf a Single Warhead 



pattern 

circle 


Assume that the target area is homogeneous and Infinitely large 
and that a single warhead bursts above target area at location (4, il). The 
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location of a target element inside the pattern circle of radius R Is 

described in polar coordinates by (r, 6) centered on (^, ti). Define 
to be the probability of kill of the target element, condi¬ 
tioned on the target element’s location as described above. Next, define fn 
to be the pdf of the rocket ballistic dispersion, and therefore il) is the 

probability that the warhead bursts at location (4, t^). Now define 
f*^[(r,0)(^,77)] to be the unconditional probability of kill of a target element 

whose location is given by (r,e) as follows; 

Next, define PktC^.'H) to be the unconditional probability of kill by a single 
warhead which bursts at location (4, "nl and is given by: 

2icJ? 

0 0 

2itK 

0 0 

Finally, define the Mean Area of Effectiveness (MAE) as: 

2nK 

MAE = J J Pj„r[(r,e)/{lrj)]rdrde (4) 

0 0 

The Mean Area of Effectiveness is commonly used to compare effective¬ 
ness of an artillery warhead. Using the MAE, we have: 

PM^,Tj)=:MAEf„(^,v) (5) 

b. Cluster Bomblet Warhead—Lethal Area 

Consider the special case where the target elements are 
uniformly distributed in a circle of radius R, and a cluster of m bomblets 
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(an equivalent area element that Is killed by one bomblet). Then is 
the probability that a target element will be killed by at least one bomblet 
in the warhead and is given by: 


'KW 



( 6 ) 


c. Evaluation of the Proportion of Target Area That is 
Destroyed by a Skdvo of n Rockets 

All the dispersion factors can be divided into two categories: 

• Precision factors that create a dispersion of the salvo around the 
Mean Point of Impact (MPI) 

• Accuracy factors that create a bias miss distance error of the vdiole 
salvo from the target center 

The bias and precision factors are independent and ran¬ 
domly distributed according to the following: 

Precision: f^iX,Y) ~ BIVARIATE NORMAL (o.O, ) 

Bias: t]) ~ BIVARIATE NORMAL (0,0, <7^.<t,) 

Using the following diagram. 


Mean Impact Point 

noget Area Unit 


target 












we can define the following probabilities: 

= Probability of killing a target area element at (a,P) with one 
warhead. 
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MAE 


e 


(7) 







and 

= Probability of killing a target area element at (a.p) with a salvo 
of n warheads. 




( 8 ) 


Recall from equation 6 that is the probability that at least 
one bomblet from a warhead of m bomblets will kill the target. TTie term 
Pg in equation 8 will use PJ^ in the Target Damage Assessment Program 
given in Appendix B. Additionally, we define as the expected 

beaten area of the target covered after a vdiole salvo: 

4, ii]. c(i, d)=^ 


This is illustrated in the following diagram: 



Generalizing, we have: 


£[-v]=rr^ 




'd^dr] 


( 10 ) 
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where: 


E[Aj\ The expected kill area from all possible salvos with all possible 

Xo.Yo The center of all possible salvos that can be taken as the cen¬ 

ter of the target (approximately, because sometimes a small 
fixed bias exists). 

For a salvo of n rounds, we have: 



=iLtl 

A-a-b 


( 11 ) 


where 

The proportion of target that is killed by a salvo of n rounds, 
d. Numerical Solution Method 

The above relationships an be evaluated in the following 


manner: 



A large space around the target, of size 3Gjd> x 3oyb. is 
divided into small cells. The center of each cell will be denoted as the 
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argument (^."n) and the corresponding G(^,ti) will be calculated for each 
cell. Then, 

s[‘^]=ZSc(«.n)/.(«.)?)AJr,Ay. (12) 

«ci Iicr 


where 






fe] ] 


(13) 


and the proportion of the destroyed target or the probability that one ele¬ 
ment target Is destroyed by a salvo of n rounds is: 


=£bl 


^ab 


(14) 


2. Model Application 

A FORTRAN program was written Incorporating the above 
model. The FORTRAN program code Is given In Appendix A The following 
results are from one set of system data made up for this program. 
Results: 

DESTRUCTION PROBABILITY = .4036 
DESTROYED AREA (SQM) = 2018. 

PDF INTEGRAL = .9946 

Data: 

NO. OF ROUNDS = 10 

MEAN AREA OF EFFECTIVENESS (SQM) = 7877.5 
TARGET HALF WIDTH (M) = 100.00 
TARGET HALF LENGTH (M) = 50.00 
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PRECISION CROSS-RANGE STANDARD DEVIATION (M) = 100.0 
PRECISION RANGE STANDARD DEVIATION (M) = 100.0 
BIAS RANGE STANDARD DEVIATION (M) = 100.0 
BIAS CROSS RANGE STANDARD DEVIATION (M) = 100.0 
NO. OF SEGMENTS ON TARGET 10. X 10. 

NO. OF SEGMENTS ON THE OPEN AREA 25.X 25. 

The target pattern that Is used for this program Is a quarter of a 
rectangle. Since the target Is homogeneous, the killing probability of each 
area element In one quarter Is Identical to a corresponding element In the 
other three. Therefore, It Is sufficient to Integrate over one-quarter of a 
target rather than the entire rectangle. 

3. Program Execution Analysis 

1. The number of segments in the open area (last line) were found to 
be insensitive— any values greater than 10 gave very similar results. 

2. PDF Integral (line 3) is an Integration of equation 13 over the whole 
space. It Is given as an indication of the accuracy of this numerical 
integration. The value accepted is sufficiently close to 1. 

3. The above results will be used in comparison with the Monte Carlo 
simulation that will be described later. 

B. SIMULATIONS OF TARGET DAMAGE ASSESSBftENT FROM 
CLUSTER BOMBLET WARHEADS 

1. Simtilation Program of Small Targets When Firing Unit Size 
is Not Considered 

The program reads system and target data. It calls for the Salvo- 
round subroutine to simulate the damage level caused by firing variable 
rounds in a salvo. 
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While receiving each one of the salvo results, the program calls 
for the Statistics subroutine that calculates averages, standard devia¬ 
tions. and confidence Intervals consecutively. 

The program stops when the confidence Interval for 95 percent 
confidence level Is reduced to predetermined fixed values (“resolution"). 
Simulation of five target types that are deployed on the same area can be 
handled simultaneously. 

Typically, any combination of five target types can be chosen 
from the following: tank. APC. artillery gun, soft vehicle, standing troops, 
prone troops, and troops in excavations. 

a. Subroutine Salvoround 

The subroutine accepts the target elements linked by a 
linked list matrix, generates bias error, generates precision error, 
searches for target elements inside the pattern ellipse— if foimd, they are 
eliminated from the list— and counts the number of eliminations. When 
the number of predetermined rounds in the salvo are completed, the 
subroutine returns the number of destroyed target elements to the main 
program. 

The program uses subroutine Rannum. which is a random 
number generator. (Professor M. Bail^ of the Naval Postgraduate School 
distributed this program to students of his Simulation class; the 
FORTRAN code list is given in Appendix B.) 

b. Assumptions 

1. The bomblet pattern on the ground is an ellipse in which the bomb- 
lets are dispersed evenly. 

2. Bias distribution and precision distribution are both independent 
bivariate normal with all means equal to 0. 
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c. Testing the Program 

Running the program using the made-up data shown in 
subsection A. 2 of this chapter with the same target yielded the following 
result: ‘‘damage level" (called “destruction probability" or “destroyed area 
proportion") » 0.397. Comparing it to 0.4036. which is obtained by the 
integration program, there is a difference of 1.6 percent. That difference 
can be explained by the termination of the simulation program when the 
half 95 percent confidence interval reaches a predetermined value, which 
had been previously fixed at 1.6 percent. Therefore, the results of the two 
programs agree (with a 95 percent confidence level). 

Figures 1 through 3 describe the procedure. 

The FORTRAN program code is given in Appendix B. 

2. Simulation Program for Small Targets When Firing Unit Size 
and Rate of Firing Are Considered 

This program is similar to the previous one. The difference is 
that the bomblet lethal area, when firing on troop target elements, is 
decreased during the elapsed salvo time. 

It is conceivable that standing troops, while exposed to fire, will 
try to find shelter. At first, they go to a prone position, and then they 
crawl to a shelter, which is assumed to be an excavation. For the purpose 
of the simulation, it Is assumed that all troops are standing at time zero 
(i.e., the lethal area values are those of standing troops) and the lethal 
area is decreased exponentially with elapsed salvo time. 

The lethal area at time t is given by: 

ALT = ALX + (ALST - ALX) • EXP (- GAMMA • TIME) 
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READ TARGET DATA 



Figure 1. Simulation Alain Program Plow Chart 
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Figure 3. Subroutine Statistics Flow Chart 

where GAMMA is the Individual-finding-shelter rate and ALST and ALX 
are lethal areas of troops when standing and when in excavations, 
respectively. 

C. SDfULATION OP TARGET DABXAGE ASSBSSblBNT FROM SMART 
SUBMUNITION WARHEADS 

This simulation assessed damage to targets consisting of hard ele¬ 
ments (armor vehicles, tanks) when hit by Sensor Fuzed Munitions 
(SFM). 

The main program is similar to the main program for cluster muni¬ 
tions (Bl); it calls for subroutine Smartscdvo. This subroutine accepts the 
target elements (listed in x. y location arrays), generates a bias error, 
generates precision error, generates submunition dispersion, and then 
searches for target elements inside the submunition searching circle. If 
targets are found, they are recorded in a stack for detected elements. 
Among the potentially detected targets, one killed target per one 
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submunition Is chosen. In order of the farthest from the searching center 
first. Finally, the subroutine counts the number of killed targets. When 
the predetermined number of salvo rounds Is completed, the subroutine 
returns the number of destroyed units to the main program. A flow chart 
of the subroutine Is given In Figure 4. 

The following two additional computations are required: 

1. Function PHI: 

This function determines the hit probability, given the existence of a 
target Inside the submunition searching circle. 

2. Function PKl: 

This function determines the kill probability, given the existence of 
a target inside the submunition searching circle. 

The FORTRAN program codes are listed in Appendix B. 

D. OTHER PROGRAMS 

1. Dealing With Large Targets 

When the target is large. It Is necessary to shift the aiming point 
In order to cover the entire target area. 

The multi-alming point procedure is essential, particularly with 
a more accurate weapon. A more accurate weapon hits the target at the 
center and the margins remain uncovered. If the weapon ballistic disper¬ 
sion is high, this procedure is less Important because the dispersion 
creates target coverage. 

Functions Splltx and Splity deviate the aiming point for large 
targets from the center to the outer portions of the target. These func¬ 
tions evenly allocate the salvo rounds among the subtargets. Residuals 
are directed to central portions. 
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Figure 4. Subroutine Sxxuurtsalvo Plow Chart 
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FORTRAN function codes are listed In Appendix B. 

2. Consecutive Statistical Calculations 

The following function and subroutines were also used: 

• Subroutine Statistics calculates the average of a current incre¬ 
mented sample and calls for subroutine Stdev to calculate the stan¬ 
dard deviations. 

• Subroutine Stdev calculates the standard deviations. 

• Function Confint calculates the confidence intervals. 

• Subroutine Rannum: A random number generator that provides all 
basic distributed random numbers. 

FORTRAN codes are listed in Appendix B. 

E. RUNNING THE PROGRAM 
1. Systems Data 
a. Range 

The analysis was done with a range of 27 km, which was 
chosen in order to include the 155 mm extended-range howitzer gun 
projectile (whose maximum range is 27 km) in the comparison analysis. 
Despite the general opinion of “ejqjerts” that modem guns can reach 
ranges longer than 27 km, the author’s opinion is that this range is the 
practical limit of standard artillery. Tube artillery will remain at this 
range because any optional extension of range will be accomp£uiled by 
degradation of accuracy as well as rate of fire. Hence, free-fllght artillery 
rocket systems will do the job more effectively. 

The maximum range of the MARS is 33 km. The range 
could actually be longer, but it was limited to 33 km due to the high dis¬ 
persion of the rockets. 
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The TCAR range Is 45 km. The range can be extended by 
using larger caliber rockets. The accuracy of the TCAR is not dependent 
on the range. 

The main reasons this thesis used a 27 km range are: (1) 
this range is common to the above three systems; (2) by convention, 
27 km is regarded as the typical aiUlleiy range for analysis purposes; 
and (3) the range is coincident with the limit of battlefield real-time intel¬ 
ligence (that limit may be extended in the near future, when the use of 
drones (RPVs) will be more extensive). 

The following data are related to the 27 km range, 
b. Accuracy 

Bias (sometimes called “ballistic") distribution (i.e., shift of 
the whole salvo from the center of target) and precision distribution (i.e., 
dispersion of the salvo around the center of Impact) are bivariate normal 
with means 0 and standard deviations Obx. oby for bias and opx, opy for 
precision. 

The standard deviations were calculated while taking into 
account all dispersion factors. The ballistic behavior of a free-fllght object 
(a projectile as well as a rocket) creates independently distributed bias 
and precision deflections in the X and Y directions. These are shown in 
Table 1. 

The analysis was done by running the simulation program 
for cluster warheads on the targets to be described. 
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TABLE 1 


DISPERSION STANDARD DEVIATIONS 


System/Dispers. 

opxm 

OpYni 



G155 

170 

36 

108 

34 

MAFIS 

214 

187 

127 

107 

TCAR 

60 

60 

5 

5 


The target location errors are assumed to have a bivariate 


normal distribution with mean 0 and standard deviations ox and oy, 
chosen to be: 0, 25, 50 and 100 m. The Target Location Error (TLE) was 
pooled with the bias deviation as follows: 

^BXTor ~ ^^TLEX^^Brrorr ~ 

Note: The variances involved are from independent random variables. 

c. Rate of Fire 

• One gun—G155 can fire at the rate of three rounds per minute. 
Maximum rate of a battery of six— 18 rounds per minute. 

• MARS— one target is engaged by no more then one launcher at a 
time, with rate of 36 rounds per minute. 

• TCAR— 36 rockets are available on a launcher, but the ground con¬ 
trol system is capable of controlling 20 rockets in flight simultane¬ 
ously (no matter from where they were fired). Since the flight 
trajectory time Is approximately one minute, the rate of fire is 20 
rounds per minute. 

d. Marginal Price 

Cost emsilysls comparisons of target damage levels between 
the three systems were carried out. This was done by taking into consid¬ 
eration the marginal cost— the cost of the rounds that were allocated to 
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the particular target (because they directly cause the damage to the 
target). 

Another aspect of this topic is that the costs of Self-Pro¬ 
pelled Guns (SPGs) and Multiple-Rocket Launchers (MRLs) are almost 
the same. Ground support equipment and logistics are also identical. 

The exact costs of all system parts were not available. In 
addition, costs of weapon systems are subject to changes according to 
various contract factors. Therefore, reasonable costs were estimated 
based on the author's experience. Table 2 shows cost estimates for one 
round of each system. 


TABLE 2 

ROUND COSTS IN DOLLARS 


System/Cost 

Cluster Bomblet Wh. 

■cJTTiSHi 

”""^G155"™'. 

$2,000 

$37,000 

MARS 

$7,000 

$112,000 

TCAR 

$11,000 

$116,000 


•It Is assumed that SFM target production cost Is $35,000 and the canis¬ 
ters that contain It cost the same as the bomblets In the cluster version. 


e. Cluster Warhead Data 

The three systems with cluster warheads accommodated 
with CL2130 at/am bomblets are IMI Ltd. products. The bomblet Is 
equipped with a time-operated self-destruction mechanism that elimi¬ 
nates the problem of duds. 

Diameter: 42 mm 
Length: 55 mm 
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Weight: 300 grams 

Maximum hard steel penetration depth: 110 mm 

Remark: This bomblet is not the M42-M77 commonly used by the 
US and NATO. 


TABLES 

SYSTEMS CLUSTER WARHEAD DATA 


System 

Diameter (mm) 

# of Bomblets 

Radius of Pattern (m) 

G155 

155 

49 

75 

MARS 

160 

104 

100 

TCAR 

160 

104 

100 


Bomblet lethal area: 

Standing troops: 197 sqm 

Prone troops: 83 sqm 

Troops in an open excavation: 1 sqm 

(It Is assumed that only a direct hit Into the excavation is effective.) 


Tank (T62): 
APC (BTR50): 
SPG: 


3.5 sqm 
5.0 sqm 

7.0 sqm (in this case, neutralization is considered 
rather than destruction) 


Trucks: 20 sqm 

f. SFM Submunition Warhead 

Warheads are accommodated with 155 mm SFM, and each 


is carried inside a canister. Data for SFM submunitions are classified, so 


this thesis will not specifically identify the submunitions. Performances 
that are listed are the best guess of the author. 
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TABLE 4 


SYSTEMS SFM WARHEAD DATA 


System 

Diameter (mm) 

# of Canisters 

Ejection Distance (m) | 

G155 

155 

1 • " 

mean 0 

sd. 5 

MARS 

160 

3 

mean 40 

sd. 10 

TCAR 

160 

3 

mean 40 

sd. 10 


There are three stages to the operation, with three different 
probabilities. Given a target inside the searching circle: 

• Detection probability Is expected to be high, therefore considered 
to be 1. 

• Hit probability (PH)—Given the target is detected, the probability 
that the forged accelerated penetrator hits the target in a vulnerable 
location. 

• Kill probability (PK)— Given a target is hit, the probability of pene« 
tration and destruction. 

PH and PK are functions of the distance of the target from the center of 
the searching circle (see Table 5): 


TABLE 5 

SFM SUBMUNITION DATA 


Section 

Radius (m) 

PH 

PK-TANK 

PH-APC 1 

Inner Circle 

25 

0.9 

0.3 

CO 

d 

Intermediate Circle 

45 

0.65 

0.25 

0.50 

Max. Search Circle 

75 

0.50 

0.2 

0.40 


False targets were not taken into account because informa¬ 


tion about that phenomenon was not available. 
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2. Simulation Output 

The programs were run on the mainframe computer (IBM 3033), 
controlled by exec programs. The number of rounds In the salvo were 
incremented after each simulation inside the main program, going from 
minimal to maximal relevant salvo size. The output file is a matrix in 
which the first row contains the salvo size and the other rows contain the 
simulation results. The results were analyzed by a GRAFSTAT software 
package and are presented In the next chapter. 
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m. TARGET DAMAGE ASSESSMENT 

A. TARGET DESCRIPTIONS AND RESULTS 

The tools that were developed in the last chapter were implemented 

to assess damage to some typical targets. 

1. Target Number 1: Uniform And Homogeneous Density- 
Troops 

* * 

* * 

* * 

* * 

* * 

i * 

* * 

4t 

* * 

* * 

Figure 5. Target No.l 



a. Description: 

Object: Personnel 

Number of units: 100 

Target area size (m): 200(width)x200(depth) 


Positions: 

1.1 

1.2 


Standing 

Prone 






1.3 Position change during salvo: Initially troops are 

standing; then all are prone or entering into top 
open shelters. “Prone function" Is an exponen¬ 
tially deca)dng function (rate = 60 seconds). 

Target location error (one standard deviation): 

Intelligence quality: 

Good: 25 (m) 

Regular: 50 (m) 

Fair: 100 (m) 

Ammunition type: Cluster bomblet warhead (or projectile) 

b. Results of 1.1 and 1,2 

The simulation results for target number 1 standing troops 
and marginal costs are shown In Figures 6 and 7, with TLE = 0, 25, 50, 
100. The simulation results for target number 1 prone troops and 
marglnsil cost are shown In Figures 8 and 9, with TLE = 0, 25, 50, 100. 

The following results can be observed from the figures: 

1. The more accurate the weapon, the more sensitive It Is to TLE. 

2. When high target damage Is required, sensitivity to TLE Increases. 

3. Tube artillery Is usually less expensive, but the rate of kill Is 
extremely low (divide the number of rounds by the rate of fire— three 
rounds per minute). In practice, more than one gun should fire at 
the same target. 

4. MARS tends to be relatively more effective when the TLE Is large 
(more than 100 m). However, when the TLE Is higher than 100, all 
results are poor, 

5. TCAR Is much superior to the other two systems, having a high rate 
of kill and a capability of achieving a high damage level, but with a 
slightly higher cost than the G155. 
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Figure 6. Target No. 1. Standing Troops 





















damage level 
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Figure 8. Target No. 1. Prone Troops 










damage level 
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Figure 9. Target No. 1, Prone Troops 










c. Simulation of a Troop Target That Changes Position 

During Salvo 

Troops find shelter In trenches during a salvo (the rate of 
finding shelters is 60 seconds). The lethal area of one bomblet when 
troops are in shelters is one square meter. Figure 10 shows the simula¬ 
tion results. 

This case illustrates the degradation of the effectiveness of 
tube artillery when it is necessary to produce massive fire power on a 
target. The massive fire power is essential because the objects will proba¬ 
bly not stay in place. 

These results show that one rocket launcher is capable of 
producing the same amount of target damage as six guns (maximum 
49 percent damage). When firing units on the battlefield are scarce, artil¬ 
lery rocket systems are advantageous, but the cost is high. 

Again, the TCAR is far superior to the other systems (maxi¬ 
mum 74 percent after 10 rounds). 
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Figure 10. Target No. 1, Standing Troops Finding Shelters 












2. Target Number 2: Troops in an Open Area, Nonsymmetrical 
Homogeneous Target 



600 M 

Figure 11. Target No. 2 

a. Description: 

Object: Personnel 

Number of units: 75 

Target area size (m): 600 (width) x 200 (depth) 

Position: Standing, prone 

Target location error (one standard deviation): 

Intelligence quality: 

Regular: 50 (m) 

Ammunition type: Cluster bomblet warhead (or projectile) 

b. Results 

Teirget number 2 ;s a “wide" area target, so the target was 
divided into sections and the aiming point for each was determined. This 
split in aiming was done the function SPLIT, which allocates a number 
of rounds to each section proportionally. The figures were rounded to 
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integers by putting a “floor” on the outer sections and a "celling" on the 
Inner sections. 

The simulation results of target number 2 for standing 
troops and marginal costs are shown In Figure 12, with TLE = 0, 50. Fig¬ 
ure 13 shows the results for prone troops. 

Aiming points deviations in x direction are: 

155 How. gun (m): 75, -75, 225, -225 

LARI60 and TCAR (m): 0, 200, -200 

The following results can be observed from the figures: 

1. Tube artillery Is the most cost-effective, but rate of kill is poor. 

2. 'Tie artillery rocket is the least cost-effective, but the rate of kill is 
nigh— it needed the same number of rounds as the gun, but the 
rate of fire is 12 times faster. 

3. TCAR has a much higher rate of kill. It is less cost-effective than 
guns, but a much higher damage level can be achieved. 

4. Results are less sensitive to 'TLE (intuitively it is explained by the 
width of the target). 
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Figure 12. Target No. 2, Standing Troops 






damage level, 7LE - 0 MARGINAL COST PER DAMAGE LEVEL 



Figure 13. Target No. 2, Prone Troops 







3. Target Number 3: Troops in a Company at a Defense Site 



Figure 14. Target No. 3 

a. Description 

Object: Personnel 

Number of units: 76 

Target area size (m): 100 (width) x 100 (depth) 

Position: Troops in open excavations 

Target location error (one standard deviation): 

Intelligence quality: 

Regular: 50 (m) 

Ammunition type: Cluster bomblet warhead (or projectile) 

Figure 15 shows the results. 

b. Results 

The following results can be observed from the figures: 

1. The relations between system performances are the same as before. 
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Figure 15. Target No. 3, Troops in Open Shelters 














110M 


2. All -'Us are very poor—at most six percent damage. Conse¬ 
quently. the cluster bomblet as a submunition is not effective 
against fortified targets. 

We might get better results if we used a greater number of 

smaller bomblets (M77) in the warheads. 

4. Target Number 4: Self-Propelled Armored Howitzer Gun Bat¬ 
tery in Firing Deployment 



^ S«l-Prap«l«dOun(SPO) 

* SP Q A m munllon and Prepallani 
Q Batttfy Aimnunilon Dump 
. SPQCrtwMambar 
O Command Sta Paiaonnal 


Figure 16. Target No. 4 
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a. Description 


Objects: 


Number of Units: 

1.1 

Armored self-propelled guns 

6 

1.2 

Guns ammunition and propellant piles 

6 

1.3 

Crew members (standing around the gun) 

5x6 

1.4 

Battery command site (crouching personnel) 

4 

1.5 

Battery ammunition dump in cases 

2 

Target area size (m): 200 (width) x 100 (depth) 


SPG: 

“W"-shaped deployment with command post at 
the center 


Ammunition: Beside each SPG, two battery dumps 

Target location error (one standard deviation): 

Intelligence quality: 

Regular 50 (m) 

Ammunition type: Cluster bomblet warhead (or projectile) 

The following assumptions are made: 

1. Only a direct bomblet hit on ammunition or propellant, in cases or 
in the open, causes detonation. Therefore, the vulnerable area of the 
ammunition dumps is taken as their base area. 

For gun ammunition and propellant piles— 16 sqm. 

For battery ammunition dump in cases— 36 sqm. 

2. The lethal area of SPG when hit by a bomblet is low. However, it is 
more practical to use the neutralization area instead. It is sufficient 
to break down one subsystem inside the SPG in order to stop it 
from functioning. This area was taken as 10 sqm. 

3. A crew is considered to be standing troops. Command site person¬ 
nel are considered to be crouching troops and the lethal area is 
assumed to be 100 sqm. 
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b. Results 


Simulation results on gun crew members and command 
site are shown in Figure 17. Results for the SPGs, gun ammunition, and 
propellant piles and battery ammunition dump are shown in Figure 18. 
Marginal costs for the last three are given in Figure 19. 

The teirget is “wide", so two aiming points were chosen: 

Deviations in X direction are: +50 m, -50 m 

The following results can be observed from the figures: 

1. Results for the crews of SPG and command sites show similar 
behaviors as before. It is important to note that only a few TCAR 
rounds are needed to destroy the command post crew, resulting in 
neutralization of the whole battery. 

2. TCAR is superior to the other two systems, with a very high killing 
rate (see Figure 18, all three objects). 

3. Guns are the most cost-effective, but they achieve only a limited 
destruction probability due to their low rate of fire. 

4. The probability of detonating ammunition is significant. Such an 
event causes fatal results and significantly damages guns and 
crews. Therefore, it is app8irent that a cluster weapon is effective as 
a counterfire weapon. If delivered by TCAR it is extremely effective. 
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Battery command site destruction proportion of crew casualties 



Figure 17. Target No. 4, Artillery Gun Battery. TLB * 50 
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Figure 18. Target No. 4, ArtUleiy Gun Battery. TLE « 50 
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5. Target Number 5: Soft Vehicles, Logistic Convoy 



Figure 20. Target No. 5 

a. Description 

Object: Soft vehicles, trucks carrying troops 

Number of units: 10, each carrying 8 troops 

Target area size (m): 600 (width) x 25 (depth) 

Position: On road, heading X direction 

Target location error (one standard deviation): 

Standing line and good intelligence 25 (m) 

On the move and good intelligence 250 (m) along the road 

25 (m) across the road 

Passengers are sitting in the trucks, two in the cabin, six in two rows— 
at the cargo platforms, on two parallel benches. 

Ammunition type: Cluster bomblet warhead (or projectile) 

Supplementary data: 

1. Trucks lethal area— 20 sqm 

2. Passengers, in sitting position— 100 sqm. They are partially pro¬ 
tected by the truck, but they are sitting Wgh above ground level. 

b. Results 

Figure 21 shows the simulation results and marginal cost 
of damage to trucks and personnel, when the convoy is stationary. 
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VEHICLES damage level. 7LE = 25 MARGINAL COST PER DAMAGE LEVEL 



Figure 21. Target No. 5, Vehicle Convoy 












Figure 22 shows the simulation results when the convoy is 
on the move. The two right graphs show the progression of damage level 
over time. 


The following results can be observed from the figures: 

1. When the target is stationaiy, we get results similar to those from 
previous targets. The highest killing rate is attributed to the TCAR; 
moreover, the G155 is the most cost-effective. 

2. When the target is on the move, the rate of kill becomes a dominant 
factor. A comparison has been made with a firing battery of G155 
(as well as a single gun). The right two graphs of Figure 18 show 
that all three systems can produce almost the same damage pro¬ 
gression through time for both trucks and personnel. The G155 bat¬ 
tery was found to be the most cost-effective. This depends on the 
availability of enough firing units in a busy arena to engage the tar¬ 
get. If this is not the case, then the G155 is not effective at all. 
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6. Target Niunber 6: Blized-Type Target In Defense Position 
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Figure 23. Target No. 6 


a. Description 


Objects: 

Tanks 

APCs 

Personnel 


Number of Units: 
3 
3 

60 


Target area size (m): 100 m radius, round fortified site 

Position: Tanks and APCs in firing trenches; personnel 

are prone behind a sand battery surrounding the 
site. 
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Target location error (one standard deviation): 

Intelligence quality: 

Regular 50 m fThls Is a typically stationary target which 

Is usually detected properly.) 

Ammunition type: Cluster bomblet warhead (or projectile) 

SFM submunition warhead (or projectile) 


Supplementary data: 

llie lethal area for a troop lying behind a sand battery Is approxi¬ 
mately half that of prone In the open— 40 sqm. 

b. Results 

Figure 24 shows simulation results and marginal costs of 
the damage level on tanks. Figure 25 shows the same on APCs. In both 
Agures, the upper two graphs show results of simulations of the three 
systems that contain SFM submunition warheads. The lower graphs 
show the same with cluster bomblet warheads. 

The following results can be observed from the figures: 

1. Rates of kill with SFM submunitions (upper-left graphs) are high for 
the three systems. TCAR is far superior to the others. MARS and 
G155 are similar In damage per number of rounds, which means 
the rate of kill of the MARS Is much higher than that of the G155 
(because the firing rate of the MARS is much higher). 

2. Ftates of kill with cluster bomblet warheads (lower-left graphs) show 
low results for TCAR £ind poor results for the others. 

3. Cost-effectiveness with SFM submunitions (upper-right graphs): 
TCAR Is the most cost-effective, significantly higher than the 
G155. MARS yields very poor results. 

4. There is not much difference In cost-effectiveness between the dus¬ 
ter warheads and the SFM warheads, but the comparison is made 
at a low level of damage (the ordinates scales represent number of 
destroyed tanks/APCs). Note that the cluster weapon cannot 
achieve the same high levels of damage as the SFM within a rea¬ 
sonable time. 
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25. Target No. 6, APCs. SPM vs. Cluster Warheads 













5. SFM submunitions: the cost to kill the first tank with TCAR is 
$770,000 ($348,000-APC): with G155 It is $1,250,000 ($666,000- 
APC) and with MARS is $3,584,000 ($2,016,000-APC). 

Remarks: 

a. The above is not represented In the graphs. 

b. The first tank to be destroyed, on average, cannot be achieved 
with a cluster bomblet weapon within a reasonable number of 
rounds. 

6. The cost to kill the second tank with TCAR Is $1,318,000; doing so 
with MARS or with G155 is not reasonably achievable. 

7. The cost of killing the first armored vehicle (tank/APC) with TCAR Is 
$232,000. This Is a better measure of cost because the effectiveness 
of the system against tanks depends on the existence of APCs in the 
target area. If an SFM submunition attacked an APC, it would not 
be available to attack a tank. In terms of effectiveness against tanks 
only, this particular submunition is wasteful. 

The first armored vehicle for G155 would cost $407,000 and for 
MARS would be $1,120,000. 

Troops: The results for prone troops In the target yield similar rela¬ 
tionships between the systems as before, so results are not shown 
on the graphs. 
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7. Target Number 7: Armor in an Offensive Site 



Figure 26. Target No. 7, Armor in an Offensive Site 
a. Description 

Objects: Number of Units: 

Tanks 5 

APC 3 

Target area size (m): 360 (width) x 160 (depth) 

Position: Tanks and APCs in firing sites 

Target location error (one standard deviation): 

Intelligence quality 
Regular 50 m 
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Ammunition type; SFM submunition warhead (or projectile) 

Aimpoints were shifted to achieve the best 
results for each system. 

b. Results 

Figure 27 shows simulation results and marginal costs of 
the damage level on tanks and APCs. 

The following results can be observed from the figures: 

1. The two upper graphs of Figure 27 are for tanks. The lower graphs 
are for APCs. The difference in the line curvatures can be explained 
as a result of having only three APC units vs. five tanks. The mar¬ 
ginal effectiveness and cost decrease as the remaining number of 
live targets becomes small. 

2. TCAR has a far higher advantage over the other two systems in kill 
rate. G155 has an extremely low kill rate. 

3. TCAR has the best cost-effectiveness for tanks and equal cost-effec¬ 
tiveness compared to G155 for APCs. 

4. The cost to kill the first tank with TCAR is $928,000 ($694,000— 
APC). with G155 is $1,665,000 (592,000-APC), and with MARS is 
$4,032,000 ($1,792,000-APC). 

5. The cost to kill the first armored vehicle (tank/APC) with TCAR is 
$348,000, with G155 $407,000. and with MARS $1,232,000. 

8. Target Number 8: Artillery Weapon with 8PM Submunition 
as an Anti-Tank Weapon 

The target is a single tank positioned in the center of the target 
area. This is an attempt to test the SFM submunition when delivered to a 
target by means of an artillery carrier and used as a pinpoint accurate 
anti-tank weapon. The question is, can tha .weapon replace the anti-tank 
missile? 
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Figure 27. Target No. 7, SFM Warhead 











Looking at the left graphs of Figure 28, we can conclude that, of 
the three systems we are considering, the TCAR is the only feasible sys¬ 
tem. It can achieve a 50 percent kill probability with 9 to 12 rockets with 
TLE of 25 to 50, respectively. However, the costs to achieve a 50 percent 
kill probability are between $1,000,000 and $1,300,000. Figure 29 shows 
the marginal cost of TCAR for TLE 25 and 50 m, in larger scale. 

The simple answer to our question is, definitely, the SFM can¬ 
not replace anti-tank missiles because it more expensive and less rapid. 

It must be understood that when a target contains many 
objects, the stochastic process is in play. Tliat Is, in a densely populated 
target area, if an artillery weapon misses object a, it is probable that it 
will instead hit object b, thereby accomplishing Its mission (the probabil¬ 
ity rising with the object/density of the target). 
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Figure 28. Target No. 8, SFM as an Anti-Tank Weapon 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


The TCAR, as a version of a modem artillery rocket system, has all 
the artillery rocket system’s advantages. However, by enhancing accu¬ 
racy, a breakthrough in artillery concepts has been created. It is no 
longer a purely statistical weapon—it both is accurate and has massive 
kill capabilities that did not exist before. 


A. EFFECTIVENESS AND MARGINAL COST 

With a cluster bomblet warhead: 

• TCAR Is capable of performing the highest rate of kill, far above the 
others, at a relatively affordable cost as compared to tube artillery. 

• Tube artillery is the most cost-effective among the three systems. 
The TCAR is a close second, and MARS is the least cost-effective (far 
below). 

• A battery of 155 mm hoAvltzer guns will accomplish all missions 
better then a free-flight rocket system. 

With an SFM warhead: 

• The TCAR Is the only system suitable for this kind of submunition. 
One must consider the difficulties of accommodating “smart" sub¬ 
munitions into projectiles. The severe launching shocks can damage 
the submunition and moderate performances (rate of kill/range) 
cause the G155 to be unsuitable for carrying “smart" submunitions. 
The end results are not worth the technological difficulties required 
to develop this type of weapon. Furthermore, the MARS capabilities 
are poor. In this regeird, MARS carrying SFM is entirely infeasible. 


B. ACCURACY 

Since the TCAR’s accuracy does not depend on the range (to a limit). 
It has a longer maximum range then any other known artillery weapon. It 
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can hit targets located deep in an opponent’s territory or engage a wider 
frontier section. 

C. RECOMMENDATIONS 

Based on the analysis of this thesis, the author recommends that 
IMI Ltd. enter into full-scale development of the TCAR. The system is 
unique, so it has a high likelihood of being acquired. 

This thesis provides a set of programs which can be used to evaluate 
damage assessment on more sets of targets with a variety of artillery 
tem versions and updated real data. The tools developed here can be 
used to generate input data for some well-known and frequently used 
combat analysis programs. 

The analysis shows a high sensitivity of results to TLE. When it 
becomes large, the effectiveness decreased very rapidly. It must be 
emphasized that all the impressive qualities of the TCAR represented in 
this thesis depend on the availability of a real-time, accurate intelligence 
system— one that provides the firing unit with information about enemy 
locations and one that observes the target continuously, giving reports In 
real time about the tsirget’s condition and the firing results. Without 
high-quality and rapid intelligence, the potential accuracy and lethality of 
the TCAR is wasted. To accomplish this, it would be necessary to inte¬ 
grate into the firing unit some type of a Remotely Piloted Vehicle (RPV) to 
be launched by, controlled by, and report directly to the firing unit. 
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APPENDIX A 

NUMERICAL INTEGRATION OP THE PROBABILITY DENSITY 
FUNCTION OVER AN HOMOGENEOUS TARGET AREA- 

PROGRAM CODE 


PROGRAM MAIN 

C- 

C Assessment of damage level an homogeneous target subjected to 
C Artillery fire of cluster bomblet projectiles or artillery 
C rockets. 

C The program prints, as output, the probability of a target 
C (area) element to be destroyed by a salvo of n rockets. 

C- 

C INDICES: Y - RANG, X - CROSS RANGE, T - TARGET, B - BIAS 
CP- PRECISION 

C- 

C PKl-The probability to kill a target element by a single round 

C PKN-The probability to kill a target element by a salvo of N 

C round 

C XT, YT - Target coordinates 
C XB, YB - Salvo bias coordinates 
C AX, AY - Halves of rectangle target dimensions 
C DXT, DYT - Area elements of the target 

C DXB, DYB - Area elements of the bias coordinates set 

C MAE - Mean area of effectiveness. 

C NDXT, NDYT - Target number of segments 
C NDXB, NDYB - Bias coordinates number of segments 
C SIGMXP, SIGMYP - Precision, standard deviation of the rounds 
C SIGMXB, SIGMYB - Accuracy, standard deviation of the salvo 
C bias. 

C N - Number of rounds in the salvo 

C$INCLUDE: 'SYSTEM.DEF' 

CSINCLUDE: 'TARGET.DEF' 

C @0@ 

INCLUDE 'SYSTEM DEF' 

INCLUDE 'TARGET DEF' 

OPEN ( UNIT-1, FILE-'SYSTEM.DTA', STATUS-'OLD') 

OPEN ( UNIT-2, FILE-'TARGET.DTA', STATUS-'OLD') 

OPEN ( UNIT-3, FILE-'TARGET.OUT', STATUS-'OLD') 

@00 

OPEN ( UNIT-1, FILE-'/SYSTEM DATA', STATUS-'OLD') 

OPEN ( UNIT-2, FILE-'/TARGET DATA', STATUS-'OLD') 

C 

READ (1 , *) N, SIGMXP, SIGMYP, SIGMXB, SIGMYB, MAE 
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READ (2 , *) AX, AY, NDXT, NDYT, NDXB, NDYB 

CALL SUBTARGET 

STOP 

END 

SUBROUTINE SUBTARGET 
C 

C- 

REAL DP, DB, DXT, DYT, DXB, DYB, EAT, CDF, XB, YB, XT, YT 
REAL PXl, PKN, G(30,30), GA(30,30), F(30,30), PKTN(10,10) 
INTEGER I, J, K, L 
C$INCLUDE: 'SYSTEM.DEF' 

CSINCLUDE: 'TARGET.DEF' 

C @00 

INCLUDE 'SYSTEM DEF' 

INCLUDE 'TARGET DEF' 

OPEN ( UNIT-3, FILE-'/TARGET OUT’, STATUS-'OLD') 

C 

DP = 2*3.14159*SIGMXP*SIGMyP 
DB - 2*3.14159*SIGMXB*SIGMYB 
DXT - 2*AX / NDXT 
DYT = 2*AY / NDYT 
DXB - 3*SIGMXB / NDXB 
DYB - 3*SIGMYB / NDYB 
C 

EAT - 0. 

CDF - 0. 

C 

DO 100 I-1,NDXB 
DO 200 J-1,NDYB 
GA(I,J) - 0. 

XB - DXB/2. +(I-1)*DXB 
YB - DYB/2. +(J-1)*DYB 
DO 300 K-1, NDXT 
DO 400 L-1, NDYT 

XT - DXT/2. +(K-1)*DXT - AX 
YT - DYT/2. +(L-1)*DYT - AY 
PXl-(MAE/DP)*EXP((((XB-XT)/SIGMXP)**2 
& +((YB-YT)/SIGMYP)**2)/(-2.)) 

PKTN(K,L) - 1-(1-PX1)**N 
GA(I,J) - GA(I,J) + PKTN(K,L) 

400 CONTINUE 

300 CONTINUE 

G(I,J) - GA(I, J)*DXT*DYT 

F(I,J) - (l./DB)*DXB*DYB*EXP(((XB/SIGMXB)**2 
& +(YB/SIGMYB)**2)/(-2.)) 

EAT - EAT + G(I,J)*F(I,J) 

CDF - CDF + F(I,J) 

200 CONTINUE 
100 CONTINUE 
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PKN - EAT / (AX*Ay) 

WRITE (3,10) PKN, EAT, 4*CDF, N ,MAE, AX, AY 
WRITE (3,20) SIGMXP, SIGMYP, SIGMXB, SIGMYB 
WRITE (3,30) NDXT, NDYT, NDXB, NDYB 
WRITE (*,10) PKN, EAT, 4*CDF, N ,MAE, AX, AY 
WRITE (*,20) SIGMXP, SIGMYP, SIGMXB, SIGMYB 
WRITE (*,30) NDXT, NDYT, NDXB, NDYB 
10 FORMAT (' DESTRUCTION PROBABILITY- *,F6.4 , 

& /• DESTROYED AREA (SQM)-',F8.0 ,/ ' PDF INTEGRAL-',F6.4 , 
& /’ NO. OF ROUNDS-',14 ,/ ' MEAN AREA OF EFFECTIVENESS', 

& ' (SQM)-', F8.1 , /' TARGET HALF WIDTH (M)-',F7.2 , 

& /' TARGET HALF LENGTH (M)-',F7.2 ) 

20 FORMAT(' PRECISION CR STANDARD DEVIATION (M)-',F5.1 , 

& /' PRECISION RANGE STANDARD DEVIATION (M)-',F5.1 , 

& / ' BIAS RANGE STANDARD DEVIATION (M)-',F5.1 , 

& / ' BIAS CRROS RANGE STANDARD DEVIATION (M)-' ,F5.1 ) 

30 FORMAT ( ' NO. OF SIGMENTS ON TARGET *,F4.0,'X',F4.0,/ 

& ' NO. OF SIGMENTS ON THE OPEN AREA ',F4.0,'X',F4.0 /) 

RETURN 
END 



SYSTEM DEF 


REAL SIGMXP, SIGMYP, SIGMXB, SIGMYB, MAE 
INTEGER N 

COMMON/SYSTEM/N, SIGMXP, SIGMYP, SIGMXB, SIGMYB, MAE 
TARGET DEF 

REAL AX, AY, NDXB, NDYB, NDXT, NDYT 
COMMON/TARGETDEF/ AX, AY, NDXB, NDYB, NDXT, NDYT 
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APPENDIX B 


TARGET DAMAGE ASSESSBIENT. 
SIMULATION PROGRAMS CODE 


PROGRAM MAIN300 

C SIMULATION OF TARGET DAMAGE ASSESSMENT, CLUSTER BOMBLET NAREAHEADS 

C THE PROGRAM READS SYSTEM AND TARGET DATA. IT CALLS FOR SALVOROND 
C SUBROUTINE TO SIMULATE THE DAMAGE LEVEL CAUSED BY FIRING VARIABLE 
C ROUNDS IN A SALVO. 

C WHILE RECEIVING THE ONE SALVO RESULTS, THE PROGRAM CALLS FOR 
C STATISTICS SUBROUTINE THAT CALCULATES CONSECUTIVELY, AVERAGE, 

C STANDARD DEVIATIONS AND CONFIDENCE INTERVALS. 

C THE PROGRAM STOPS WHEN THE CONFIDENCE INTERVAL FOR 95% CONFIDENCE 
C IS REDUCED TO A FIXED RESOLUTION. 

C SIMULATION OF 5 TARGET TYPES THAT ARE DEPLOYED IS THE SAME AREA 
C CAN BE HANDLED SIMULTANEOUSLY. 

C TYPICALLY ANY COMBINATION OF 5 CAN BE CHOSEN OUT OF THE FOLLOWING: 

C TANK, APC, ARTILLERY GUN, SOFT VEHICLE, STANDING TROOPS, PRONE 
C TROOPS, LAYING TROOPS AND CROUCHING (IN FOXHOLES ) TROOPS. 

C 

C AL(I)...LETHAL AREA OF A BOMBLET ON TARGET TYPE I 
C AP...BOMBLET PATTERN ELLIPSOID, CROSS RANGE RADIUS 
C BP...BOMBLET PATTERN ELLIPSOID, RANGE RADIUS 
C NBOM...NUMBER OF BOMBLETS IN ONE ROUND. 

C TD(I)...NUMBER OF TARGET TYPE I ELEMENTS THAT WHERE DESTROYED. 

C RESOL...RESOLUTION 
C NEXTd, J) . . .LINK LIST MATRIX. 

C 

C- 

C$INCLUDE: 'CTARGET.DEF' 

C$INCLUDE: 'CSEED.DEF' 

C$INCLUDE: 'CSYSTEM.DEF' 

C e@0 

INCLUDE 'CTARGET DEF' 

INCLUDE 'CSEED DEF' 

INCLUDE 'CSYSTEM DEF' 

C 

INTEGER I, J, M, COUNT, TD(5), NBOM, NEXT(5,100) 

INTEGER MINSALVO, MAXSALVO, INTERVAL 

REAL TDAVG(5), TDSTD(5), RESOL(5), CONFINT, AL(5), DMGLEVEL(5,3) 
REAL TLEX, TLEY 
LOGICAL DONE 
CHARACTER*8 TTYPE(5) 

C 

C 003 

OPEN ( 1, FILE-'/TARGETS DATA', STATUS-'OLD') 

OPEN ( 2, FILE-'/WEAPON DATA', STATUS-'OLD') 

OPEN ( 3, FILE-'/MAIN300 OUT', STATUS-'OLD') 
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OPEN ( 4, FILE-'/CSEED DATA', STATUS-'OLD' ) 
OPEN ( 5, FILE-'/MAIN300 LIST', STATUS-'OLD') 


READ ( 2, *) SIGXP, SIGYP, SIGXB, SIGYB, RATE, AP, BP, NBOM 
READ ( 2, *) ( AL{I), I-l, 5 ) 

READ (2, *) NPLAT, MINSALVO, MAXSALVO, INTERVAL, LTIME 
READ ( 2, *) NAPX, NAPY, DX, DY, TLEX, TLEY 
READ ( 4, *) { SEED(I), 1-1,6) 

DATA RESOL / .1, .1, .1, .1, .6 / 

SIGXB - SORT(SIGXB**2 + TLEX**2) 

SIGYB - SQRT(SIGYB**2 + TLEY**2) 

DO 100 I-l, 5 

READ (1, *) N(I) , TTYPE(I), ( DMGLEVEL{I,M), M - 1, 3 ) 
DO 110 J-1, N(I) 

READ (1, *) X(I,J), Y{1,J) 

NEXT(I, J) - J + 1 
110 CONTINUE 

IE ( N(I) .GT. 0 ) NEXTd, N(I)) - 0 

PKl(I) - 1 - ( 1 - AL(I) / ( 3.14159*AP*BP ) )**NBOM 
100 CONTINUE 

WRITS (5, 40) SIGXB, SIGYB, NAPX, DX, TLEX, TLEY 

DO 200 NSALVO - MINSALVO, MAXSALVO, INTERVAL 
COUNT - 0 
DO 201 I-l, 5 
TD(I) - 0 
TDAVG(I) - 0.0 
201 CONTINUE 


300 CONTINUE 

PRINT*, 'FROM MAIN COUNT-',COUNT, 'AP-',AP, 'BP-',BP 
DONE - .TRUE. 

CALL SALVOROUND ( NEXT, NSALVO, TD ) 

DO 350 I - 1, 5 

IF ( N(I) .GT. 0 ) THEN 

CALL STATISTICS (COUNT, TD(I), TDAVG(I), TDSTD(I) ) 
IF(CONFINT(COUNT,TDSTD(I)).GT.RESOL(I))DONE -.FALSE. 
IF ( COUNT .LT. 3 ) DONE - .FALSE. 

END IF 

350 CONTINUE 

COUNT - COUNT + 1 
IF ( .NOT. DONE ) GO TO 300 


WRITE (3, *) 

WRITE (3, *) ' SALVO-',NSALVO , ' REPEAT-',COUNT 
WRITE (3 ,10) ( TTYPE(I), TDAVG(I), I - 1, 5 ) 

WRITE (3 ,20) ( TTYPE(I), TDAVG(I)/(N(I)+.001), 1 - 1, 5) 
WRITE (5 ,30) NSALVO,((TDAVG(I)/(N(I)+.001))*100, I - 1, 5) 
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200 CONTINUE 

C WRITE (*, *) 'PKKD- ',( PKl(I), 1-1,5 ),' 
C 
C 
C 

10 FORMAT (/,5(1X,A7,•-•,F4.1) ) 

20 FORMAT { 5(1X,A7,'-•,F4.3) ) 

30 FORMAT ( IX,14,5(IX,F4.1) ) 

40 FORMAT (' SIGXB-',F7.2,' SIGXB-',F7.2,' NAPX' 
& • TLEX-',F6.2, ' TLEy-',F6.2 ) 

STOP 

END 


REPEAT-',COUNT 


',12,' DX-',F6.1, 
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SUBROUTINE SALVOROUND( NEXTMATRIX, NSALVO, SCORE ) 

C SIMULATION OF DAMAGE LEVEL WHEN FIRING A SALVO OF FIXED ROUNDS 
C 

C THE SUBROUTINE ACCEPTS THE TARGET ELEMENTS LINKED BY A LINK LIST 
C MATRIX. SAMPLE A BIAS ERROR. SAMPLE PRECISION ERROR. SEARCH FOR 
C A TARGET ELEMENT IN THE PATTERN ELLIPSOID. IF FOUND, IT IS 
C ELIMINATED FROM THE LIST. COUNT NUMBER OF ELIMINATIONS. 

C WHEN SALVO PREDETERMINED ROUNDS WHERE COMPLETED, RETURNS. 

C 

C- 

C 

CINCLUDE: 'CTARGET.DEF' 

CINCLUDE: 'CSEED.DEF' 

CINCLUDE: 'SYSIOO.DEF' 

C @00 

INCLUDE 'CSYSTEM DEF' 

INCLUDE 'CTARGET DEF' 

INCLUDE 'CSEED DEF' 

C @00 

INTEGER ROUND, I, J, SCORE(5), FIRST(5), PNT, NEXT(5,100) 

INTEGER NEXTMATRIX(5, 100), NSALVO 

REAL TIME, XB, XH, XP, YB, YH, YP , RAN, SPLITX, SPLITY 
LOGICAL FST 
C 

CALL RANNUM( 2, SEED(l), 0.0, SIGXB, 0, XB ) 

CALL RANNUM( 2, SEED (2), 0.0, SIGYB, 0, YB ) 

C 

ROUND - 0 
TIME - 0.0 
DO 101 I - 1, 5 
FIRST(I) - 1 
SCORE(I) - 0 
DO 102 J - 1, N(I) 

NEXT(I, J) - NEXTMATRIX(I, J) 

102 CONTINUE 
101 CONTINUE 
C 

100 CONTINUE 

ROUND - ROUND + 1 

CALL RANNUM( 2, SEED (3), 0.0, SIGXP, 0, XP ) 

CALL RANNUM( 2, SEED (4), 0.0, SIGYP, 0, YP ) 

XH - XB + XP 

YH - YB + YP 

IF ( NAPX .GT.l ) XH - XH + SPLITX(ROUND, NSALVO, NAPX, DX) 

IF ( NAPY .GT.l ) YH - YH + SPLITY(ROUND, NSALVO, NAPY, DY) 

DO 300 I - 1, 5 

IF ( N(I) .GT. 0 ) THEN 
150 CONTINUE 

IF ( FIRST(I) .GT. 0 ) THEN 
FST - .FALSE. 

J - FIRST(I) 

C WRITE(3,*) 'FROM SALVO 1 ROUND-',ROUND,’ J-',J 

IF((((XH-X(I, J) )/AP) ** 2 + ((YH-Y(I, J) )/BP) **2) .LT. 1 .)THEN 
CALL RANNUM( 1, SEED (5), 0.0, 1.0, 0, RAN ) 
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IF ( RAN .LT. PKl(I) ) THEN 
FIRST(I) - NEXT(I,J) 

SCORE(I) - SCORE(I) + 1 
FST - .TRUE. 

END IF 
ENDIF 

IF ( FST ) GO TO 150 
PNT “ FIRST(I) 


200 

200 IF ( NEXT(I, J) .GT. 0 ) THEN 

J - NEXT(I, J) 

IF ( ( ( (XH-X (I, J) ) /AP) **2+ ( (YH-y (I, J) ) /BP) **2) .LT. 1 .) THEN 
CALL RANNUM( 1, SEED (5), 0.0, 1.0, 0, RAN ) 

IF ( RAN .LT. PKl(I) ) THEN 
NEXT(I,PNT) - NEXT(I,J) 

SCORE (I) «■ SCORE (I) + 1 
ELSE 

PNT - J 
ENDIF 
ELSE 

PNT - J 
ENDIF 
GO TO 200 
ENDIF 
ENDIF 
ENDIF 

300 CONTINUE 


IF ( ROUND .LT. NSALVO ) GO TO 100 

RETURN 

END 


INCLUDE FILES 
CSYSTEM DEF 
C 

INTEGER NAPX, NAPY 

REAL SIGXP, SIGYP, SIGXB, SIGYB, RATE, AP, BP, PK1{5), DX, DY 
COMMON /SYSTEM/ SIGXP, SIGYP, SIGXB, SIGYB, RATE, AP, BP, PKl, 

& NAPX, NAPY, DX, DY 
C 

C PKl(I)...PROBABILITY TO KILL A TARGET TYPE I IN THE PATTERN OF ONE 
C WARHEAD 

C RATE...RATE OF FIRE OF ONE PLATFORM 
C NAPX, NAPY... NUMBER OF AIMPOINTS ON THE TARGET 

C DX, DY ... MAX. DISTANCE OF AIMPOINT DEVIATION FROM THE TARGET CENTER 
C 


CTARGET DEF 
C 

INTEGER N(5), NPLAT 
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c 


REAL X(5,100), Y(5,100), LTIME 
COMMON/TARGET/ X, Y, N, NPLAT, LTIME 


CSEED DEF 
C 

INTEGER SEED(8) 
COMMON/SEED SET/ SEED 
C 
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PROGRAM MAINTIME 


C 

C SIMULATION OF TARGET DAMAGE ASSESSMENT, CLUSTER BOMBLET WAREAHEADS 
C IN THIS SIMULATION LETHAL AREA DECREASES AS A FUNCTION OF SALVO 
C TIME DURATION. 

C --— 

c 

C THE PROGRAM READS SYSTEM AND TARGET DATA. IT CALLS FOR SALVOROND 
C SUBROUTINE TO SIMULATE THE DAMAGE LEVEL CAUSED BY FIRING VARIABLE 
C ROUNDS IN A SALVO. 

C WHILE RECEIVING THE ONE SALVO RESULTS, THE PROGRAM CALLS FOR 
C STATISTICS SUBROUTINE THAT CALCULATE CONSECUTIVELY, AVERAGE, 

C STANDARD DEVIATIONS AND CONFIDENCE INTERVALS. 

C THE PROGRAM STOPS WHEN THE CONFIDENCE INTERVAL FOR 95% CONFIDENCE 
C IS REDUCED TO A FIXED RESOLUTION. 

C SIMULATION OF 5 TARGET TYPES THAT ARE DEPLOYED IS THE SAME AREA 
C CAN BE HANDLED SIMULTANEOUSLY. 

C TYPICALLY ANY COMBINATION OF 5 CAN BE CHOSEN OUT OF THE FOLLOWING: 
C TANK, APC, ARTILLERY GUN, SOFT VEHICLE, STANDING TROOPS, PRONE 
C TROOPS, LAYING TROOPS AND CROUCHING (IN FOXHOLES ) TROOPS. 

C 

C AL(I)...LETHAL AREA OF A BOMBLET ON TARGET TYPE 1 
C AP...BOMBLET PATTERN ELLIPSOID, CROSS RANGE RADIUS 
C BP...BOMBLET PATTERN ELLIPSOID, RANGE RADIUS 
C NBOM...NUMBER OF BOMBLETS IN ONE ROUND. 

C TD(I)...NUMBER OF TARGET TYPE I ELEMENTS THAT WHERE DESTROYED. 

C RESOL...RESOLUTION 
C NEXT(I,J)...LINK LIST MATRIX. 

C 

C- 

C$INCLUDE: 'CTARGET.DEF’ 

C$INCLUDE: 'CSEED.DEF' 

CSINCLUDE: 'C400SYST.DEF' 

C @00 

INCLUDE 'CTARGET DEF' 

INCLUDE 'CSEED DEF' 

INCLUDE 'C400SYST DEF' 


C 


C 

C 

C 

C 

C 

C 000 


C 


C 


INTEGER I, J, M, COUNT, TD(5), NEXT(5,100) 
INTEGER MINSALVO, MAXSALVO, INTERVAL 
REAL TDAVG(5), TDSTD(5), RESOL(5), CONFINT 
LOGICAL DONE 
CHARACTER*8 TTYPE(5) 


OPEN ( 1, 
OPEN ( 2, 
OPEN ( 3, 
OPEN ( 4, 


FILE-'TARGET.DTA', STATUS-'OLD') 
FILE-'SYSTEM.DTA', STATUS-'OLD') 
FILE-'SYSTEM.OUT', STATUS-'OLD') 
FILE-'SEED.DTA', STATUS-'OLD' ) 


OPEN 

OPEN 

OPEN 

OPEN 

OPEN 

OPEN 


( 1, FILE-'/TARGET 1 DATA', STATUS-'OLD') 
( 2, FILE-'/WEAPON DATA', STATUS-'OLD') 

( 2, FILE-'/CSYSTEM DATA', STATUS-'OLD') 
( 3, FILE-'/TIME OUT', STATUS-'OLD*) 

( 4, FILE-'/CSEED DATA', STATUS-'OLD' ) 

( 5, FILE-'/TIMELIST OUT', STATUS-'OLD') 
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READ ( 2, *) SIGXP, SIGYP, SIGXB, SIGYB, RATE, AP, BP, NBOM 
READ ( 2, *) ( AL(I), I-l, 5 ) ,AL0 

READ ( 2, *) NPLAT, MINSALVO, MAXSALVO, INTERVAL,(GAMMA(I),I-l,5) 
READ ( 4, *) ( SEED(I), 1-1,6) 

C 

DATA RESOL / .8, .8, .8, .8, .8 / 

C 

DO 100 I-i, 5 

READ (1, *) N(I) , TTYPEd) 

DO 110 J-1, N(I) 

READ (1, *) X(I,J), Y{I,J) 

NEXT(I, J) - J + 1 
110 CONTINUE 
C 

IF ( N(I) .GT. 0 ) NEXTd, Nd)) - 0 
100 CONTINUE 
C 

DO 200 NSALVO - MINSALVO, MAXSALVO, INTERVAL 
COUNT - 0 
DO 201 I-l, 5 
TDd) - 0 
TDAVGd) - 0.0 
201 CONTINUE 


300 CONTINUE 

PRINT*, 'FROM MAIN COUNT-',COUNT, ’AP-',AP, 'BP-',BP 
DONE - .TRUE. 

CALL SALVOTIME( NEXT, NSALVO, TD ) 

DO 350 I - 1, 5 

IF ( Nd) .GT. 0 ) THEN 

CALL STATISTICS ( COUNT, TDd), TDAVGd), TDSTDd) ) 
IF(CONFINT(COUNT,TDSTD(I)).GT.RESOL(1))DONE -.FALSE. 
IF ( COUNT ,LT. 3 ) DONE - .FALSE. 

END IF 

350 CONTINUE 

COXJNT - COUNT + 1 
IF ( .NOT. DONE ) <30 TO 300 


WRITE (3, *) 

WRITE (3, *) ' SALVO-',NSALVO , ' REPEAT-',COUNT 
WRITE (3 ,10) ( TTYPEd), TDAVGd), I - 1, 5 ) 

WRITE (3 ,20) ( TTYPEd), TDAVG (1) / (N (1) +. 001) , I - 1, 5) 
WRITE (5 ,30) NSALVO, ((TDAVGd)/(Nd)+.001)) *100, I - 1, 5) 

200 CONTINUE 


10 FORMAT (/,5(1X,A7,'-•,F4.1) ) 
20 FORMAT ( 5(1X,A7,'-',F4.3) ) 
30 FORMAT ( IX,14,5(IX,F4.1) ) 
STOP 
END 
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SUBROUTINE SALVOTIME( NEXTMATRIX, NSALVO, SCORE ) 


C SIMULATION OF DAMAGE LEVEL WHEN FIRING A SALVO OF FIXED ROUNDS 
C IN THIS SIMULATION LETHAL AREA IS DECREASED AS FUNCTION OF THE 
C SALVO TIME DURATION 
c 
C 

C THE SUBROUTINE ACCEPTS THE TARGET ELEMENTS LINKED BY A LINK LIST 
C MATRIX. SAMPLE A BIAS ERROR. SAMPLE PRECISION ERROR. SEARCH FOR 
C A TARGET ELEMENT IN THE PATTERN ELLIPSOID. IF FOUND, IT IS 
C ELIMINATED FROM THE LIST. COUNT NUMBER OF ELIMINATIONS. 

C WHEN SALVO PREDETERMINED ROUNDS WHERE COMPLETED, RETURNS. 

C 

C- 

c 

CINCLUDE: 'CTARGET.DEF' 

CINCLUDE: 'CSEED.DEF* 

CINCLUDE; 'C400SYST.DEF' 

C @@@ 

INCLUDE 'C400SYST DEF» 

INCLUDE 'CTARGET DEF' 

INCLUDE 'CSEED DEF' 

C @00 

INTEGER ROUND, I, J, SCORE{5), FIRST(5), PNT, NEXT(5,100) 
INTEGER NEXTMATRIX(5, 100), NSALVO 
REAL TIME, XB, XH, XP, YB, YH, YP , RAN, PKl, ALT 
LOGICAL FST 
C 

CALL RANNUM( 2, SEED(l), 0.0, SIGXB, 0, XB ) 

CALL RANNUM( 2, SEED(2), 0.0, SIGYB, 0, YB ) 

C 

ROUND - 0 
TIME - 0.0 
DO 101 I - 1, 5 
FIRST(I) » 1 
SCORE(I) " 0 
DO 102 J - 1, N(I) 

NEXT(I, J) - NEXTMATRIX(I, J) 

102 CONTINUE 
101 CONTINUE 
C 

100 CONTINUE 

ROUND = ROUND + 1 

TIME - TIME + ROUND / NPLAT / RATE 

CALL RANNUM( 2, SEED(3), 0.0, SIGXP, 0, XP ) 

CALL RANNUM( 2, SEED(4), 0.0, SIGYP, 0, YP ) 

XH - XB + XP 
YH = YB + YP 
DO 300 I - 1, 5 

IF ( N(I) .GT, 0 ) THEN 

ALT - AL0+(AL(I)-AL0)*EXP(-GAMMA(I)*TIME) 

PKl - 1. - ( 1. - ALT / (3.14159*AP*BP ) )**NBOM 
150 CONTINUE 

IF ( FIRST(I) .GT. 0 ) THEN 
FST - .FALSE. 
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J » FIRST(I) 

C WRITE (3,*) 'FROM SALVO 1 ROUND-ROUND,' J-',J 

IF((((XH-X(I,J))/AP)**2+((YH-y(I,J))/BP) ** 2 ) .LT.1.)THEN 
CALL RANNUM( 1, SEED(5), 0.0, 1.0, 0, RAN ) 

IF ( RAN .LT. PKl ) THEN 
FIRST(I) - NEXT(I,J) 

SCORE(I) = SCORE(I) + 1 
FST = .TRUE. 

END IF 
ENDIF 

IF ( FST ) GO TO 150 
PNT - FIRST(I) 

C 

C 200 
C 

200 IF ( NEXT(I, J) .GT. 0 ) THEN 

J = NEXT(I,J) 

IF{ ( ( (XH-X(I, J) ) /AP) **2+{ (YH-Yd, J) ) /BP) **2) .LT.1.)THEN 
CALL RANNUM( 1, SEED(5), 0.0, 1.0, 0, RAN ) 

IF ( RAN .LT. PKl ) THEN 
NEXT(I,PNT) = NEXT(I,J) 

SCORE(I) = SCORE(I) + 1 
ELSE 

PNT - J 
ENDIF 
ELSE 

PNT - J 
ENDIF 
GO TO 200 
ENDIF 
ENDIF 
ENDIF 

300 CONTINUE 
C 
C 

IF { ROUND .LT. NSALVO ) GO TO 100 
C 

RETURN 

END 


INCLUDE FILE 
C400SYST DEF 
C 
C 

REAL SIGXP, SIGYP, SIGXB, SIGYB, RATE, AP, BP, AL{5), GAMMA(5) 
REAL ALO 
INTEGER NBOM 

COMMON /SYSTEM/ SIGXP, SIGYP, SIGXB, SIGYB, RATE, AP, BP, 7 ,, 

& NBOM, GAMMA, ALO 
C 

C RATE...RATE OF FIRE OF ONE PLATFORM 
C 

CTARGET DEF 

INTEGER N(5), NPLAT 
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REAL X(5,100), Y(5,100), LTIME 
COMMON/TARGET/ X, Y, N, NPLAT, LTIME 
C 

CSEED DEF 
C 

INTEGER SEED(8) 

COMMON/SEED SET/ SEED 
C 
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PROGRAM SMARTMAIN 

C SIMULATION OF TARGET DAMAGE ASSESSMENT, SMART SUBMUNITION WAREAHEADS 
C 

C THE PROGRAM READS SYSTEM AND TARGET DATA. IT CALLS FOR SMARTSALVO 
C SUBROUTINE TO SIMULATE THE DAMAGE LEVEL CAUSED BY FIRING VARIABLE 
C ROUNDS IN A SALVO. 

C WHILE RECEIVING THE ONE SALVO RESULTS, THE PROGRAM CALLS FOR 
C STATISTICS SUBROUTINE THAT CALCULATE CONSECUTIVELY, AVERAGE, 

C STANDARD DEVIATIONS AND CONFIDENCE INTERVALS. 

C THE PROGRAM STOPS WHEN THE CONFIDENCE INTERVAL FOR 95% CONFIDENCE 
C IS REDUCED TO A FIXED RESOLUTION. 

C SIMULATION OF 5 TARGET TYPES THAT ARE DEPLOYED IS THE SAME AREA 
C CAN BE HANDLED SIMULTANEOUSLY. 

C TYPICALLY ANY COMBINATION OF 5 CAN BE CHOSEN LIKE THE FOLLOWING: 

C TANK, APC, SELF-PROPELLED ARTILLERY GUN, SOFT VEHICLE. 

C 

C TD(I)...NUMBER OF TARGET TYPE I ELEMENTS THAT WHERE DESTROYED. 

C RESOL...RESOLUTION 
C 

C- 

C$INCLUDE: ’CTARGET.DEF' 

C$INCLUDE: 'CSEED.DEF' 

C$INCLUDE: 'SMARTSYS.DEF' 

C$INCLUDE: 'SEARCH.DEF' 

C 

INCLUDE 'CTARGET DEF' 

INCLUDE 'CSEED DEF' 

INCLUDE 'SMARTSS DEF' 

INCLUDE 'SEARCH DEF' 

C 

INTEGER I, J, COUNT, TD(5) 

INTEGER MINSALv'O, MAXSALVO, INTERVAL, NSALVO 
REAL TDAVG(5), TDSTD(5), RESOL(5), CONFINT, TLEX, TLEY 
LOGICAL DONE 
CHARACTER*8 TTYPE{5) 

C 

C OPEN ( 1, FILE-'SMARTTAR.DTA', STATUS-'OLD') 

C OPEN ( 2, FILE-'SMARTSYS.DTA', STATUS-'NEW') 

C OPEN ( 3, FILE-'SYSTEM.OUT', STATUS-'OLD') 

C OPEN ( 4, FILE-'SEED.DTA' , STATUS-'OLD* ) 

C @00 

OPEN ( 1, FILE-'/TARGET7A DATA', STATUS-'OLD') 

OPEN ( 2, FILE-'/SMARTSS DATA', STATUS-'OLD') 

OPEN { 3, FILE-'/CSYSTEM OUT', STATUS-'OLD') 

OPEN ( 'i, FILE-'/CSEED DATA', STATUS-'OLD' ) 

OPEN ( 5, FILE-'/SMARTMN LIST', STATUS-'OLD' ) 

C 

READ (2, *) SIGXP,SIGYP,SIGXB,SIGYB,SIGMAR,DIST,RSEARCH,R1,R2 
READ { 2, *)(PRH1{I), I-l, 5), (PRH2(I), 1-1,5), (PRH3{I), 1-1,5) 
READ ( 2, *)(PRK1(I), I-l, 5), (PRK2(I), 1-1,5), (PRK3(I), 1-1,5) 
READ (2, *) NPLAT, MINSALVO, MAXSALVO, INTERVAL, NSUB 
READ (2, *) NAPX, NAPY, DX, DY, TLEX, TLEY 
READ ( 4, *) ( SEED(I), 1-1,8) 

C 
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DATA RESOL / .1, .1, .1, .1, .01 / 

SIGXB - SQRT(SIGXB**2 + TLEX**2) 

SIGYB - SQRT(SIGYB**2 + TLEY**2) 

DO 100 I-l, 5 

READ (1, *) N(I) , TTYPE(I) 

DO 100 J-1, N{I) 

READ (1, *) Y{I,J) 

100 CONTINUE 

WRITE (5, 40) SIGXB, SIGYB, NAPX, DX, TLEX, TLEY 

DO 200 NSALVO - MINSALVO, MAXSALVO, INTERVAL 
COUNT - 0 
DO 201 I-l, 5 
TD(I) - 0 
TDAVG(I) - 0.0 
201 CONTINUE 


300 CONTINUE 

PRINT*, 'FROM MAIN COUNT-',COUNT, 'AP-',AP, 'BP-',BP 
DONE - .TRUE. 

CALL SMARTSALVO < NSALVO, TD ) 

DO 350 I - 1, 5 

IF ( N(I) .GT. 0 ) THEN 

CALL STATISTICS ( COUNT, TD<I), TDAVG<I), TDSTD(I) ) 
IF(CONFINT(COUNT,TDSTD(I)).GT.RESOL{I))DONE -.FALSE. 
IF ( COUNT .LT, 10) DONE - .FALSE. 

ENDIF 

350 CONTINUE 

COUNT - COUNT + 1 
IF ( .NOT. DONE ) GO TO 300 


WRITE (3, *) 

WRITE (3, *) ' SALVO-',NSALVO , ' REPEAT-',COUNT 
WRITE (3 ,10) ( TTYPE(I), TDAVG(I), 1 - 1, 5 ) 

WRITE (3 ,20) { TTYPE(I), TDAVG(I)/(N(I) +.001), 1 - 1, 5) 
WRITE (5 ,30) NSALVO,((TDAVG{I)/{N(I) +.001))*100, I - 1, 5) 
200 CONTINUE 


10 FORMAT (/,5(1X,A7,•-',F4.1) ) 

20 FORMAT ( 5(IX,A7,•-•,F4.3) ) 

30 FORMAT ( IX,14,5(IX,F4.1) ) 

40 FORMAT (' SIGXB-', F7 . 2, ' SIGYB-' ,F7.2,' NAPX-', 12,' DX-',F6.1, 
& ' TLEX-',F6.2, ' TLEY-',F6.2 ) 

STOP 

END 
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SUBROUTINE SMARTSALVO( NSALVO, SCORE ) 

C 

c-..--...— 

C THIS SUBROUTINE SIMULATES TARGET DAMAGE LEVEL AFTER FIRING A SALVO 
C OF SMART WARHEAD ROUNDS. 

C 

C 

C THE SUBROUTINE ACCEPTS THE TARGET ELEMENTS, LISTED IN X Y LOCATION 
C ARRAYS. SAMPLES A BIAS ERROR. SAMPLES PRECISION ERROR. SAMPLE 
C SUBMUNITION DISPERSION. SEARCHES FOR A TARGET ELEMENTS INSIDE THE 
C THE SUBMUNITION SEARCHING CIRCLE. IP FOUND, IT IS RECORDED IN A 
C DETECTED LIST. AMONG THE POTENTIALLY DETECTED TARGETS, ONE KILLED 
C TARGET, PER ONE SUBMUNITION IS CHOSEN, IN ORDER OF THE FARER FROM THE 
C SEARCHING CENTER, THE FIRST TO BE SAMPLED. COUNT NUMBER OF KILLED 
C TARGETS. WHEN SALVO PREDETERMINED ROUNDS IS COMPLETED, RETURNS. 

C 

C- 

c 

C$INCLUDE: 'CTARGET.DEF' 

CSINCLUDE: 'CSEED.DEF' 

C$INCLUDE: 'SMARTSS.DEF' 

C 800 

INCLUDE 'CTARGET DEF' 

INCLUDE 'CSEED DEF' 

INCLUDE 'SMARTSS DEF' 

C 000 

INTEGER ROUND, I, J, L, K, SCORE(5), 11(500), JJ(500) 

INTEGER NSALVO, SUB(500), ISUB, PNT, IPNT, ND 
REAL XB, YB, XP, YP, THETA, RAN, XCS, YCS, RR(500), SPLITX 
REAL XX(5,100), YY(5,100), RMAX, DISTR, RANGE, PHI, PKl, SPLITY 
C 

CALL RANNUM( 2, SEED(l), 0.0, SIGXB, 0, XB ) 

CALL RANNUM( 2, SEED(2), 0.0, SIGYB, 0, YB ) 

C 

IF ( NAPX .GT.l ) XB - XB + SPLITX(ROUND, NSALVO, NAPX, DX) 

IF ( NAPY .GT.l ) YB - YB + SPLITY(ROUND, NSALVO, NAPY, DY) 

C 

ROUND - 0 
DO 101 1 - 1, 5 
SCORE(I) - 0 
DO 102 J - 1, N(I) 

XX(I, J) - X(I, J) 

YY(I, J) - Y(I, J) 

102 CONTINUE 
101 CONTINUE 
C 

100 CONTINUE 

C PRINT*, 'I AM IN THE SALVO LOOP ,ROUND-',ROUND,(N(I),I-l,5) 

ROUND - ROUND + 1 

CALL RANNUM( 2, SEED(3), 0.0, 5IGXP, 0, XP ) 

CALL i<ANNUM( 2, SEED (4), 0.0, SIGYP, 0, YP ) 

PNT - 1 
C 

DO 200 ISUB - 1, NSUB 

CALL RANNUM( 2, SEED(4), DIST, SIGMAR, 0, DISTR ) 

CALL RANNUM( 1, SEED(5), 0.0, 6.28319, 0, THETA ) 
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XCS - XB + XP + DISTR*COS(THETA) 

YCS - YB + yP + DISTR*SIN(THETA) 

PRINT*, 'XCS-',XCS,'YCS-',YCS 

DO 300 I «■ 1, 5 

IF ( N(I) .6T. 0 } THEN 
DO 400 J - 1, N(I) 

RANGE - SQRT( (XCS-XX(I,J))**2+(YCS-YY(I,J))**2 ) 
IF ( RANGE .LT. RSEARCH ) THEN 
II(PNT) - I 
JJ(PNT) - J 
RR(PNT) - RANGE 
SUB(PNT) - ISUB 
PNT - PNT + 1 
END IF 


400 

CONTINUE 


ENDIF 

300 

CONTINUE 

200 

CONTINUE 


700 

CONTINUE 


RMAX - 


DO 500 L - 1, PNT - 1 

IF ( SUB(L) .GT. 0 ) THEN 
IF ( RR(L) .GT. RMAX ) THEN 
RMAX - RR(L) 

IPNT - L 
END IF 
ENDIF 

500 CONTINUE 

080000 

WRITE(3, 555) ( II(LL),JJ(LL),RR(LL),SUB(LL) ,LL-1,PNT-1) 

555 FORMAT (14,14,F6.2,13) 

WRITE(3,*)'IPNT-',IPNT,'PNT-',PNT,'RMAX-',RMAX,'ROUND-',ROUND 

000080 

IF ( RMAX .LT. 0.1 ) GO TO 800 

CALL RANNUM( 1, SEED(6), 0.0, 1.0, 0, RAN ) 

IF ( RAN .LT. PH1(II(IPNT),RMAX) ) THEN 
ND - SUB(IPNT) 

DO 600 K - 1, (PNT - 1) 

IF ( SUB(K) .EQ. ND ) SUB(K) - 0 
600 CONTINUE 

CALL RANNUM( 1, SEED(7), 0.0, 1.0, 0, RAN ) 

IF ( RAN .LT. PKl(II(IPNT),RMAX) ) THEN 

.IN CASE OF A KILLED TARGET CAN BE HIT AGAIN 

IF(II(K).EQ.II(IPNT).AND.JJ(K).EQ.JJ(IPNT).AND. 

& XX(II(K)) .EQ. 1000000 ) ) 

.MUST BE PUT AS NEGATIVE CONDITION 

SCORE(II(IPNT)) - SCORE(II(IPNT)) + 1 
XX(II(IPNT),JJ(IPNT)) - 1000000. 
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DO 650 K - 1, (PNT - 1) 

IF(II(K) .EQ.II(IPNT) .AND.JJ(K) .EQ. JJ(IPNT) ) RR(K) — 
650 CONTINUE 

END IF 
ENDIF 

RR(IPNT) - -1. 

GO TO 700 
C 

800 CONTINUE 
C 

IF ( ROUND .LT. NSALVO ) GO TO 100 

PRINT*, (SCORE(I), I«l,5) 

RETURN 
END 


REAL FUNCTION PHI(I, R) 


THIS FUNCTION DETERMINES THE HIT PROBABILITY GIVEN EXSISTANCE 
OF TAGET INSIDE THE SUBMUNITION SEARCHING CIRCLE. 


INTEGER I 
REAL R 

SINCLUDE; 'SEARCH.DEF' 

INCLUDE 'SEARCH DEF' 

IF ( R .LT. R1 ) THEN 
PHI - PRHld) 

ELSEIF ( R .LT. R2 ) THEN 
PHI- PRH2(I) 

ELSE 

PHI - PRH3(I) 

ENDIF 

RETURN 

END 


REAL FUNCTION PK1(I, R) 


THIS FUNCTION DETERMINES THE KILL PROBABILITY GIVEN EXSISTANCE 
OF TAGET INSIDE THE SUBMUNITION SEARCHING CIRCLE. 


INTEGER I 
REAL R 

C$INCLUDE: 'SEARCH.DEF' 

INCLUDE 'SEARCH DEF' 

IF ( R .LT. R1 ) THEN 
PKl - PRKl(I) 

ELSEIF ( R .LT. R2 ) THEN 
PKl- PRK2(I) 

ELSE 

PKl - PRK3(I) 

ENDIF 

RETURN 

END 
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INCLUDE FILES 


CTARGET DEF 
C 

INTEGER N(5), NPLAT 
REAL X(5,100), y(5,100), LTIME 
COMMON/TARGET/ X, Y, N, NPLAT, LTIME 
C 

CSEED DEF 
C 

INTEGER SEED(8) 

COMMON/SEED SET/ SEED 
C 

SMARTSS DEF 
C 

REAL SIGXP, SIGYP, SIGXB, SIGYB, RATE, RSEARCH, SIGMAR, DIST 
REAL DX, DY 

INTEGER NSUB, NAPX, NAPY 

COMMON /SYSTEM/ SIGXP, SIGYP, SIGXB, SIGYB, RATE, RSEARCH, NSUB, 
& SIGMAR, DIST, NAPX, NAPY, DX,DY 


SEARCH DEF 
C 

REAL Rl, R2, PRH1(5), PRH2(5), PRH3(5), PRKl(5), PRK2(5), PRK3(3) 
COMMON /SEARCHCIRCLE/ Rl, R2, PRHl, PRH2, PRH3, PRKl, PRK2, PRK3 
C 
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REAL FUNCTION SPLITX(ROUND, NSALVO, NAPX, DX) 

C THIS FUNCTION DEVIATES THE AIMPOINT FOR LARGE TARGET 
C 

INTEGER ROUND, NSALVO, NAPX, N, M 
REAL DX 

IF (NSALVO .LT. NAPX) THEN 
SPLITX - 0. 

GO TO 200 
ENDIF 

IF( INT( REAL(NAPX)/2. ) .LT. ( ( REAL(NAPX)/2.) - 0.1 ) ) THEN 
M-1 
N-0 

100 CONTINUE 

IF (ROUND .LT. ((REAL(M)/REAL(NAPX) )*NSALVO + 0.7)) THEN 
SPLITX - 2.*DX/(NAPX-1)*N*((-1)**M) 

GO TO 200 
ENDIF 
M - M + 1 

N - INT( REAL(M) / 2. ) 

GO TO 100 

ELSE 

M-1 

N-0 

120 CONTINUE 

IF ( ROUND .LT. (( REAL(M)/REAL(NAPX) )*NSALVO + 0.6)) THEN 
SPLITX - 2.*DX/REAL(NAPX-1)*( N+.5 )*((-1)**(M+1)) 

GO TO 200 
ENDIF 
M - M + 1 

N - INT( REAL(M-1) / 2. ) 

GO TO 120 
ENDIF 

200 CONTINUE 
RETURN 
END 


REAL FUNCTION SPLITY(ROUND, NSALVO, NAPY, DY) 

Q ******************************************************************* 

C THIS FUNCTION DEVIATES THE AIMPOINT FOR DEEP TARGET 
C 

INTEGER ROUND, NSALVO, NAPY, N, M 
REAL DY 

IF (NSALVO .LT. NAPY) THEN 
SPLITX - 0. 

GO TO 200 
ENDIF 

IF( INT( REAL(NAPY)/2. ) .LT. { ( REAL(NAPY)/2.) - 0.1 ) ) THEN 
M-1 
N-0 

100 CONTINUE 

IF ( ROUND .LT. (( REAL(M)/REAL(NAPY) )*NSALVO + 0.6)) THEN 
SPLITY - 2,*DY/(NAPY-1)*N*((-1)’»*M) 

GO TO 200 
ENDIF 
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M - M + 1 

N - INT( REAL(M) 12 .) 

GO TO 100 

ELSE 

M-1 

N-0 

120 CONTINUE 

IF ( ROUND .LT. (< REAL(M)/REAL(NAPT) )*NSALVO + 0.7)) THEN 
SPLITY - 2.*DY/REAL(NAPY-1)*( N+.5 )*((-1)**(M+1)) 

GO TO 200 
ENDIF 
M - M + 1 

N - INT( REAL(M-1) 12 .) 

GO TO 120 
ENDIF 

200 CONTINUE 
RETURN 
END 
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REAL FUNCTION CONFINT(REP,STD) 

C— 

INTEGER REP 
REAL STD 

IF(REP.GT.2) THEN 

CONFINT-4.* STD/SQRT(REAL(REP-1)) 

ELSE 

CONFINT - 4. 

ENDIF 

RETURN 

END 

SUBROUTINE STATISTICS (REPEAT,NEWX,AVGX,STDX) 

C THIS SUBRUTINE CALCULATES THE AVARAGE OF A CURRENT 
C INCREMENTED SAMPLE, AND CALLS FOR SUBRUTINE STDEV 
C TO CALCULATE THE STANDARD DEVIATION. 

C~»— 

INTEGER REPEAT, NEWX 
REAL AVGX, STDX, XAVOLD 
XAVOLD-AVGX 

C WRITE(*,*) 'I AM IN STATISTCS REPEAT-',REPEAT 
XAVOLD-AVGX 

AVGX-(REPEAT*AVGX+NEWX)/REAL(REPEAT+1) 

CALL STDEV (REPEAT, NEWX, AVGX, STDX, XAVOLD) 

RETURN 

END 

Q* it ******* it * ilititit** ********** it ***** it hhli** ********* K*** 

SUBROUTINE STDEV (REPEAT, NEWX, AVGX, STDX, XAVOLD) 

^s:«taK«c»aBaB«iai«tata«aBMaaMaBasasase«BaBMaKaa:aiSMaaB«rTOWaKSsaBiBMaBaifliatwaB» 

INTEGER REPEAT, NEWX 
REAL AVGX, STDX ,XAVOLD, RP 
IF (REPEAT .EQ.O) THEN 
STDX-0.0 
RETURN 
ELSE 

RP-REAL(REPEAT) 

STDX-SQRT((RP-1.)/RP*STDX**2+XAV0LD**2+NEWX**2/RP-(RP+1.)/RP 
& *AVGX**2) 

ENDIF 

RETURN 

END 
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c_ 

SUBROUTINE RANNUM(DISTN, SEED, RPARMl, RPARM2, IPARM, X) 

C__ 

C THIS SUBROUTINE PROVIDES AN INTERFACE WITH THE LLRANDOMII 

C ROUTINES PROVIDED IN THE NONIMSL LIBRARY. THE PARAMETER 

C REQUIREMENTS AND CALLING PROCEDURES ARE AS FOLLOWS: 

C 

C DISTN - DISTRIBUTION TYPE YOU WANT TO SELECT 
C AN INTEGER BETWEEN 1 AND 7 

C SEED = THE RANDOM NUMBER SEED YOU WISH TO USE 

C RPARMl, RPARM2, AND IPARM ARE REAL AND INTEGER PARAMETERS 

C PASSED TO THE ROUTINE WITH MEANINGS WHICH VARY WITH THE TYPE OF 

C DISTRIBUTION YOU DESIRE 

C NOTE: IPARM IS CURRENTLY NOT BEING USED. 

C X - THE RETURNED RANDONM NUMBER, IT IS ALWAYS REAL 
C 

C DISTRIBUTION NUMBERS AND THE ASSOCIATED PARM DEFINITIONS 

C 

C 1—UNIFORM ON THE INTERVAL RPARMl TO RPARM2 

C 2—NORMAL WITH MEAN RPARMl AND STD DEV RPARM2 

C 3—EXPONENTIAL WITH RATE RPARMl 

C 4—COUCHY WITH A » RPARMl AND B - RPARM2 

C 5—GAMMA WITH SHAPE RPARM2 AND RATE RPARMl 

C 6—POISSON WITH RATE RPARMl 

C 7—GEOMETRIC WITH P = RPARMl 

C 

C NOTE TO NEW USERS: THIS FUNCTION NOW HAS TWO CALLS TO LLRAN 

C FOR EACH TYPE OF VARIATE, ONE COMPATABLE WITH THE MAINFRAME 

C LLRANDOMII, AND ONE COMPATABLE WITH RNDG.LIB FOR THE PC 

C THE PC VERSIONS ALL END IN 'PC’ 

C ONE IS ALWAYS COMMENTED OUT. 

C 

C VERY IMPORTANT: IF YOU ARE USING THE LLRANDOMII CALLS, MAKE 

C SURE TO COMMENT OUT THE STATEMENT 

C SEED - INT{PCSEED) 

C 

REAL RPARMl, RPARM2, X 
INTEGER DISTN, SEED, IPARM, N 
DOUBLE PRECISION PCSEED 
REAL TEMP, VARIAT(l) 

TRANSLATION FOR PCs 

PCSEED - DFLOAT{SEED) 

N - 1 


IF (DISTN.LE.O.OR.DISTN.GT.8) THEN 

WRITE(10, *) 'ILLEGAL CALL TO RANNUM, BAD DISTN' 
STOP 
END IF 


WRITE(10, *) 
WRITE(10, *) 
WRITE(* , *) 
WRITE(* , *) 
WRITE(*, *) 


'ee@003@e00eeee@e@ee0e0ee0@@e0eeeeeee00ee6' 

'FROM RANNUM' 

'DISTRIBUTION NUMBER - ', DISTN 
'SEED - ', SEED 
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c 

WRITE(* , 


'PCSEED 

- ', PCSEED 

c 

WRITE(* , 

*) 

'RPARMS 

ARE ', RPARMl, ' ', RPARM2 


GOTO (10, 

20, 

30, 40, 

50, 60, 70), DISTN 


c 

C GENERATE A UNIFORM BETWEEN RPARMl AND RPARM2 

10 CONTINUE 

IF (RPARMl - RPARM2.EQ.0) THEN 
C WRITE(10, *) 'ILLEGAL EQUAL RPARMS IN RANNUM' 

STOP 
END IF 

IF (RPARM1.GT.RPARM2) THEN 
TEMP - RPARMl 
RPARMl - RPARM2 
RPARM2 - TEMP 
ENDIF 
C 

C CALL LRNDPC(PCSEED, VARIAT, N) 

CALL LRND(SEED, VARIAT, 1, 1, 0) 

C 

VARIAT(1) - RPARMl + (RPARM2 - RPARMl) * VARIAT(1) 

GOTO 99 
C 

C GENERATE A NORMAL WITH MEAN RPARMl AND STDDEV RPARM2 
20 CONTINUE 

C CALL LNORPC(PCSEED, VARIAT, N) 

CALL LNORM(SEED, VARIAT, 1, 1, 0) 

VARIAT(1) - (VARIAT(1) * RPARM2) + RPARMl 
GOTO 99 
C 

C GENERATE AN EXPONENTIAL WITH RATE (1/MEAN) RPARMl 

30 CONTINUE 

IF (RPARMl.EQ.O) THEN 

C WRITE(10, *) 'ILLEGAL ZERO RATE IN RANNUM' 

STOP 

ENDIF 

C CALL LGAMPC(PCSEED, VARIAT, N, 1.0) 

CALL LEXPN(SEED, VARIAT, 1, 1, 0) 

VARIAT(1) - VARIAT(1) / RPARMl 
GOTO 99 
C 

C GENERATE A COUCHY WITH A - RPARMl AND B - RPARM2 

40 CONTINUE 

IF (RPARM2.LE.0) THEN 

C WRITE(10, *) 'ILLEGAL COUCHY SPREAD IN RANNUM, B - ', RPARM2 

STOP 
ENDIF 

C CALL LCHYPC(PCSEED, VARIAT, N) 

CALL LCCHY(SEED, VARIAT, 1, 1, 0) 

VARIAT(1) - (VARIAT(1) * RPARM2) + RPARMl 
GOTO 99 

50 CONTINUE 

IF (RPARMl.LE.O) THEN 

C WRITE(10, *) 'ILLEGAL NONPOSITIVE GAMMA RATE IN RANNUM' 
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STOP 
END IF 

IF (RPARM2.LE.0) THEN 

• C WRITE(10, *) 'ILLEGAL SHAPE PARAMETER IN RANNUM' 

STOP 
END IF 

‘ C CALL LGAMPC(PCSEED, VARIAT, N, RPARM2) 

CALL LGAMACSEED, VARIAT, 1, 1, 0, RPARM2) 

VARIAT (1) - VARIAT (1) * (1.0 / RPARMl) 

GOTO 99 

60 CONTINUE 

IF (RPARMl.LE.O) THEN 

C WRITE(10, *) 'ILLEGAL POISSON RATE IN RANNUM' 

STOP 

ENDIF 

C CALL LPOIPC(PCSEED, VARIAT, N, RPARMl) 

CALL LPOIS(SEED, VARIAT, 1, 1, 0, RPARMl) 

GOTO 99 

70 CONTINUE 

IF (RPARMl.LE.O) THEN 
WRITE(* , *) 'RETURNING ', X 

wRiTEdo, *) 'e@eeeeeeeeeeeeseeseeeeeeeeeeee6eeeeeeee0@' 

WRITE(10, *) 

RETURN 
END 
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